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Abstract 



The gravitational instability of inhomogeneities in the expanding universe 
is studied by the relativistic second-order approximation. Using the tetrad for- 
malism we consider irrotational dust universes and get equations very similar 
to those given in the Lagrangian perturbation theory in Newtonian cosmology. 
Neglecting the cosmological constant and assuming a flat background model 
we give the solutions of the nonlinear dynamics of cosmological perturbations. 
We present the complete second-order solutions, which extend and improve 
earlier works. 
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I. INTRODUCTION 



Gravitational instability and structure formation in the universe is an important topic of 
cosmological research. By using N-Body codes it is possible to follow the general nonlinear 
evolution of initially small perturbations numerically, but an understanding of what has 
happened between input and output can often better be gained by analytical treatments. 
Various analytical approaches have been compared with the numerical results statistically 
||T|J^ and it has turned out that the celebrated Zel'dovich approximation gives the best 
fit to the numerical treatment. Buchert presented the Lagrangian perturbative approxi- 
mation to first-order based on Newtonian theory. This work was extended to second-order 
IP and even to third-order giving some new useful information about self-gravitating 
systems. But these and most other analytical treatments are Newtonian approaches, which 
are valid only for perturbations on scales much smaller than the horizon size. On super- 
horizon scales instead one needs a relativisic approach. The pioneer was Lifshitz with his 
linearized theory on the basis of general relativity, which was extended to the second-order 
by Tomita ||8|-[T0|. A Lagrangian relativistic approximation to the second-order based on 
fluid flow equations was given by Matarrese, Pantano and Saez |[TT|Jl^ . Parry, Salopek and 
Stewart |T3| presented the nonlinear solution of the Hamilton- Jacobi equation for general 
relativity, using the spatial gradient expansion technique W^ and reproduced the Zel'dovich 
approximation. The "Higher-order Zel'dovich approximation" is discussed in Croudace et 



al. and Salopek, Stewart and Croudace |l6 



In this paper, we give an alternative approach, which extends a tetrad-based Zel'dovich- 



type approximation by Kasai We derive the fully general relativistic equations very 

similar to those given in the Newtonian case, which are solved in a fiat background model 
without the cosmological constant by an iteration method. The complete solutions are 
compared with previous work and it is found that they include all these results. 

This paper is organized as follows. In Sec. II we present the basic relativistic equations 
and introduce the tetrad formalism. In Sec. HI the perturbative approach is presented and 
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the solutions up to second-order are given. In Sec. IV we compare our results with previous 
works. Sec. V contains conclusions. In the appendices we explain our gauge condition and 
present the complete second-order solutions, including the decaying and the coupling mode. 
Units are chosen so that c = 1. Indices /x, z/, ■ ■ ■ and a,b, - ■ ■ run from to 3 and indices 
■ ■ ■ run from 1 to 3. 



II. EXPOSITION OF THE METHOD 

In this section, we summarize a general relativistic treatment to describe the non-linear 
evolution of an inhomogeneous universe ||T7|-[19| . The models we consider contain irrotational 
dust with density p and four- velocity u'^ (and possibly a cosmological constant A). Neglecting 
the fluid pressure and the vorticity is a reasonable assumption in a cosmological context. In 
comoving synchronous coordinates, the line element can be written in the form 

ds^ = -dt^ + gijdx'dx^ (2.1) 

with = (1,0,0,0). Then the Einstein equations read 

i [ + {k\)' - K^^K\ ] = 8nGp + A , (2.2) 
i^^ll, - = , (2.3) 
k'^ + K\K'^ + = {AnGp + A) 6'^ , (2.4) 

where ^R^ j is the three-dimensional Ricci-tensor, 

K', = lg%k (2.5) 

is the extrinsic curvature, || denotes the covariant derivative with respect to the three metric 
gij, and an overdot (') denotes d/dt. The energy equation u^T'^^.^ = gives 

p + pK\ = 0, (2.6) 

with the solution 
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/det[^,,(t,„.,x)] 

p = p{tin, x) ' , , =r ■ (2-7) 

/det [gij{t,x)\ 



The evolution equation for the Ricci curvature is obtained in the form [|T7 

^i?^. + 2K\ 'r) = + k) \ - - II, • (2.8) 

Let us introduce the scale factor function a{t) which satisfies the following equation 

2ad + a^ + k-Aa'^ = , (2.9) 

where the curvature constant k takes the value of +1, 0, —1 for closed, flat, and open spaces, 
respectively. (Eq. ( p.9|) is obtained from the Friedmann equation 

a\2 k 8nG A 
aj 6 6 



t) +TI = ^Pfc + T (2-10) 



and its derivative with respect to time.) If the spacetime is exactly Friedmann-Lemaitre- 
Robert son- Walker (FLRW), then we have 

K'. = -5'^, ^Ry = 2^5'^ for FLRW. (2.11) 

di (X 

Therefore, the deviations from the FLRW models due to inhomogeneity are expressed by 
the peculiar part of the extrinsic curvature 

V) = K'. - -5'^ , (2.12) 

which represents the deviation from the uniform Hubble expansion, and the deviation of the 
curvature tensor 

n'j=a^^R'^-2kS'^. (2.13) 



Using these quantities, Eqs. (p.3|), ( |2.4| ), and ( p.8|) are rewritten as 

V) + (3^ + n) V) + 1 (7^^. - 5v) = \ {{vir - vWi] 5^, , (2.15) 



^ e\j ^ ^ j \e ^ j \i ^ e \j 



] ' t J ' • J • e\i ' • i \e • i \e • i n ' (2-16) 
where | denotes the covariant derivative with respect to the conformally transformed three 



metric 'jij = a gij. 



The procedure essential to develop the relativistic Zel'dovich-type approximation is 
to introduce the following orthonormal tetrad 



with 



j(0) 



= (-1,0,0,0), e(2 = (0,e(?)^(0,a(t)e?) for £ 



1,2,3. 



The spatial basis vectors are parallelly transported along each fluid line, i.e.. 



In our choice of the tetrad components, it reads 



or V\ 



Using Eqs. ( |2.15| ), ( |2.16| ) and (|2.20|) , we obtain the following key equation 



d_ 

dt 



where 



dt\ A 



yk yn 

' n. k 



- a' {V\V\ - V\V\ , e 



Q 



and 



S 



(^) 



(2.17) 

(2.18) 

(2.19) 

(2.20) 

(2.21) 

(2.22) 
(2.23) 



(2.24) 



Note that the left-hand-side of Eq. ( |2.21| ) is already linearized with respect to e^P, and 
all terms on the right- hand- side, except S'^f , are manifestly nonlinear quantities. It has, 
therefore, a form suitable for solving it perturbatively by iteration. It should also be stressed 
that we have not used any approximation methods in deriving Eq. ( ^^.21] ). Our treatment 
here is fully nonlinear and exact. 



III. PERTURBATIVE APPROACH 



In this section, we solve perturbatively the key equation Eq. (|2.21 ) by an iteration 
method. 

A. The background 

The background ( V^j = 0, Tl'j = ) solution is characterized by 

e^?=0, i.e., e(? = e(?(x). (3.1) 

Furthermore, the metric = S(^k){e)^^'i^G^^j is that of a constant curvature space with curva- 
ture constant k. In the case of a fiat background, we can write 

e^? = 5^? for k = 0. (3.2) 

Hereafter, we restrict our consideration to the Einstein-de Sitter background, /c = 0, A = 0. 

B. The first-order solutions: scalar modes 

Linear perturbations are classified into scalar, vector, and tensor modes. In the first- 
order level, they do not couple with each other, and can be discussed separately. Let us first 
consider the scalar perturbations. The general form for the linearly perturbed triad in this 
case is 

e? = 5? + i^? = 5? + 5^?(^^^^ + ^lO • (3-3) 

Let us write the first-order quantities with subscript (1). Then the perturbed extrinsic 
curvature is 

V^\),=i5V + 5'V (3.4) 
From the constraint equation (|2.14|), which reads 
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(3.5) 



in the first-order, we obtain A j = 0. However, the part A{t) 5*^- in the extrinsic curvature 
simply represents the uniform and isotropic Hubble expansion. Therefore, by a suitable 
re-definition of the background, we can set 



A = , i.e., A = A(x) 



(3.6) 



As was noted previously, it is apparent that the source terms P^P and Q^^ are second- 
order quantities (and higher). Using V^-^-jj = B'^j, we also find that S^f vanishes in linear 
order: 



C(^) _ ( yj ,k I yk ,j _ yj ,k _yk ,j \ _ n 



(3.7) 



Therefore, to first-order, the right-hand-side of the key equation (|2.21|) vanishes and it can 
be integrated to give 



(3.8) 



By choosing C^?(x) = -AnGpba^ S^f (A(x) 5{ + C%x)), Eq. (jSj) is re-written as 



dt 



(3.9) 



Note that now it has the same form as the equation which governs the density contrast 6 in 
conventional linear perturbation theory ||20]. Using the growing mode D'^{t) = a{t) = t^/^ 
and the decaying mode solutions D~{t) = respectively, we obtain the solutions in the 
form 



5'V = C^(x) + t3 ^'^^^.(x) + $'V(x) . 



(3.10) 



For the metric, we have the following first-order expression: 



Qi, = a\t) [(1 + 2A(x)) S^j + 2C,i,(x) + 2ti ^,i,(x) + 2t-^ <l>,i,(x) 



(3.11) 



The relation between A{x.) and \E'(x) is given by Eq. (|2.15|) . To first-order, it reads 



where 



^(1). + 3^^(\), + ^ fei). - ) = ' (3-12) 



2 1 - -2/F>,i 



^(\), = 3t-^n-^ (3-13) 



and 



n\,,^ = -A'\^-A'\5%. (3.14) 



''(l)i J ,k" j 

Hence we have 

A(x) = Hvl;(x) . (3.15) 

The function C(x) is not determined by the Einstein equations within our approximation. 
As shown in APPENDIX A, however, we can set C(x) = using a residual gauge freedom. 
The final form of the first-order solutions is, therefore, 

= (l + y vl;(x)) 5? + 5^'^ (t§ ^'\,{^) + $'^,(x)) , (3.16) 

or in the form of the metric 

20 



9i, = a\t) [l + y^(x)) 5ij + 2t§ ^,,,(x) + 21'^ <l>,i,(x) 



(3.17) 



Note that we have not assumed that the density contrast is small, in order to derive the 
solutions. The density is given by Eq. ([Z^ ), which in this case reads 



/, h{Un)\ det e^i^(ti„,x) 

p = p{tin, X — — — . 3.18 

V / det e7 t,x 



C. The first-order solutions: tensor modes 

Under the assumption of vanishing vorticity, the remaining is the tensor mode perturba- 
tions. In this case, we can write the triad in the form 



^ = 5^ + ^H^, (3.19) 



with H\ , = and H\ = 0. 



The perturbed extrinsic curvature is 



(3.20) 



Then, the constraint equation ( |2.14| ) to first-order, i.e., Eq. (|3.5D, is trivially satisfied. 

To obtain the equation for H^^, we can use the key equation (|2.21|) . On the right-hand- 
side of Eq. ( |2.21|) , S^f is the only quantity to be calculated since P^P and Q^? are of higher 
order: 



'^{1)1 — ° j ^ i ,k ■ 



(3.21) 



Eq. ( g:2lD reads 



d_ 

di 



Tji ,k 



(3.22) 



Integrating Eq. ( |3.22| ), we obtain 



H\ + 3-H\ - \ \/^H\ = 



(3.23) 



where is the Laplacian of fiat 3-spaces. In fact, the same equation for can be also 
obtained directly from Eq. (|3.12|) . The solution of Eq. ( p.23| ) is given as 



(fq t 2j^3(3|q|t3) /i* . exp(iq • x) 



(3.24) 



where J^s is the Bessel function of order ±3/2 and /i* ■ is a constant tensor with h^^ = 



and Qih^ j = 0. (See, e.g., Ref. for detail. 



D. The second-order solutions 

In order to avoid notational complexity, in this subsection we only deal with growing 
mode terms. The complete solutions of the decaying and coupling terms can be found in 
APPENDIX B. Moreover, we omit the first-order tensor mode. (It is not our aim to consider 
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the nonlinear effect wliicli comes from tliis mode. With respect to this problem, see Ref. 
PJTO[].) Thus we begin with the following form 



e?=(l + Hv,(x))5?+ti5(?v,.yx)+.?. (3.25) 

The second-order quantity is decomposed into a transverse-traceless part and a remaining 
longitudinal part 

= S^f {f3\ + X\) , (3.26) 

where x'j,i = 0, x\ = 0. 

The peculiar deformation tensor to second-order is immediately found to give 

Quantities with subscript (2) represent the second-order term in the expansion. From the 
constraint equation (|2.14| ), we now obtain 

/3^-.-/3\, + |H(v[/'^-vI/,,)^^. = 0. (3.28) 

Let us turn our attention to the key equation ( p.21| ). To second-order it reads 



+ 2^^'? - = lc'?(x) + 1 /' . (P,» + Qj«. + Sg!.) dt , (3.29) 



where c'-f (x) is a second-order integration "constant". It is apparent that the source terms 
P^2)i ^^'^ Q(2)j ^'^^ quadratic with respect to the first-order quantities, hence contain neither 

'(2), 



nor x^j- Furthermore, from Eq. (|3.28| ), we find that the longitudinal part of S^2)i does 



not contain Actually, if we take the divergence of Eq. (|3.29|) , we obtain 

, + 2-/3V , - 4v^Gp,/5^ , = - - ^''.^''A . . (3.30) 

J, ^ J, J, ^6 J, 3 ^ ' ' 

Therefore, solutions for can be written as a linear combination of the homogeneous 
solution and the inhomogeneous solution in the presence of the given source terms: 

= a(x) 5^ + tl V-Vx) + n <^*^.(x) , (3.31) 
10 



where we have used a convenient choice of the integration "constant", c^f(x) = 
—AnGpba^alx.) 6''j. 

Once we obtain the temporal dependency of the solutions, their spatial dependency, i.e., 
'?/'*j(x) and V^*j(x) are determined by Eq. ( p.l5| ). To second-order 

y^^- = ^ 5^- - V,) , (3.33) 

where 

^^^^ ^ \ {^>\^>^\^ - M/'^,,v^'y . (3.34) 

(The tensor /i*^ has an interesting property: the trace gives the second scalar invariant 
[] of the tensor \E'' j. ) 

The equation for j can also be obtained from Eq. ( p.29| ), but it is more convenient to 



use Eq. ( |2.15|) instead. To second-order, it gives for the transverse-traceless part 



where 



+ - = 5', , (3.35) 



'5^ = |//, + |(/.5^-V,)'S (3.36) 



is a transverse and traceless tensor: S^^ = 0, 5*^- j = 0. This shows that gravitational waves 
are induced even if there are initially scalar perturbations only. The solution of Eq. (p.35|) 



was given by Tomita |^ in the following way. Introducing the conformal time variable t], 
which is related to thy dt = adrj, Eq. ( p.35| ) is rewritten as 



^ Three scalar invariants of a three-dimensional tensor A^j are defined by I (A) = A^^, II {A) ^ 
l/2[{A\-y - A'jA\], III{A) = dei{Aj). They satisfy the relation dei{5^ ^ + A ^) = 1 + I{A) + 
II{A) + III{A). See, e.g., Ref. |p and references therein. 
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^2 



4 d 



1 



Eq. ( I3.37D can be solved using the retarded Green function as 



(3.37) 



X^(x,r/) = -j Dret(x,r/;x',V) a\i^') 5V(x') di^' <f-^ 



(3.38) 



where 



Dret (x, ; x', r/') = ^-^-^^ (1 + e(r/ - V)) (r^ - r^) + (r^ - r^)) (3.39) 



with r = |x — x'l and t = t] — t]'. Substituting Eq. ( p.39|) into Eq. (p.38|) , the solution reads 

1 rv-Vin 



X\(x,r/) 



IMAttt]^ Jo 

where r' = |x'|. 

Finally, we obtain the metric tensor up to second-order 

20 ^ 100 

+ 2a I dij 
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dr' r' {6r] + r') {r] - r') - r'r^,^ / dQ' S% (x + x') , (3.40) 



Tij ^ 9ij 

+ a{t) 



1 H ^ H + 2a (5, 



9 81 

20,2 9 \ 20 
— + -a hi + — 
9 5 9 

19_fc , 12 , fc _ .3 
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(3.41) 



We still have freedom in choosing a(x), which corresponds to the second-order term of the 
initial amplitude of the gravitational potential fluctuations. It can be absorbed into the first- 
order perturbations by a suitable re-definition of ^(x). For example, choosing a = —'^'^'^ 
gives 



+ y a(t) (-6^,,^,, - 4**,,, + m'^m^k 6,, 



(3.42) 



+ -a\t) 



+ 2x. 



At the initial time {t —* 0) only first-order metric perturbations exist. 
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IV. COMPARISON WITH PREVIOUS WORKS 

In this section, we compare our result (Eq. ( |3.42| )) with previous work. Quantities, which 
refer to these papers, will be indicated by a hat. 



A. Tomita's second-order theory 

Tomita [Q] extended Lifshitz's linearized theory |^ up to the second-order calculation on 
the basis of general relativity. Setting F = for the growing mode and F = 54<l> for 
the decaying mode (see Eq. (4.1) in Ref. his result is fully coincident with ours, except 
for one point: he did not consider the terms due to the coupling between the growing and 
decaying modes, which are included in our complete solutions. (See APPENDIX B.) 



B. The fluid flow approach 

Matarrese et al. |]12[ also carried out second-order calculations based on the fluid flow 
approach. Their result (Eq. (49) in Ref. [0]) is partly consistent with ours, since they 



neglect several terms in the computed metric. In spite of the fact that they obtain the 
initial condition from the gauge-invariant linear theory, they neglect the first-order constant 
mode, Y "^^ij our notation in Eq. ( |3.42| ) in the subsequent calculations. Also missing is 



the second-order homogeneous solution, which is proportional to t^^^. 

The comparison of the second-order transverse-traceless parts has to be taken with some 
caution. Eq. (B19) in Ref. |T^, which has to be solved, can be derived from our Eq. ( p.35|) . 



In the short-wavelength limit inside the horizon [rj 3> r') in our approach we get V^x*j = 
—t^S'^j, which can be identified with Eq. (65) in Ref. |jl2[. In the long-wavelength limit 
outside the horizon they obained a result, which can be neglected cause of the appearance 
of spatial derivatives, whereas in our exact result there exists no solution for the wavelength 
larger than the horizon size. 
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C. The gradient expansion technique 



Parry et al. |13] derived a nonlinear solution for gij, based on the gradient expansion 
method. (See also Ref. [|l^,|l^].) Their "fifth-order" result is the 



i" 9/-- -\ 81 
%j = t^kij + —t^ [Rkij - ARij) + — tS 



350 
5 



89 
32^ 



-R'^.j, — AR'^^Rhf ] hi 



5 - 

lORRij + -^R;ij + nRj^^Rjn — T^Rij-k 



(4.1) 



where kij = kij{x) is the "seed" metric, Rij and R are the 3-dimensional Ricci tensor and 
Ricci scalar, respectively, of the 3-metric kij, and a semicolon (; ) denotes the covariant 
derivative with respect to kij. To compare our solution Eq. (p.42|) with their Eq. ( [4.1|) , we 

set 



20 

= (1 + Y*W)<5^, • 



(4.2) 



Then the Ricci tensor up to second-order is 
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100^ T 



200 



vT/vl/ .. 



100 



200 



(4.3) 



If we substitute this expression into Eq. ( ^.Ij ) and calculate up to the second-order, we obtain 

20 



7ij = t 



9 

1 4 



(4.4) 



Therefore, we find that their "fifth-order" result coincides with our second-order solution, 
except for the transverse-traceless part, Xtj- If ^6 take the long- wavelength limit, Xij can be 
neglected, since spatial derivatives are assumed to be quantities of higher order than time 
derivatives in this limit and as a result, the wave equation for Xij does not appear. In this 



sense, our result includes theirs. 
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V. CONCLUDING REMARKS 



In this paper, we have developed the second-order perturbative approach to the nonhnear 
evolution of irrotational dust universes in the framework of general relativity. We have shown 
the complete calculation of the second-order solutions in a k = 0, A = background, based 
on the tetrad formalism given by Kasai [0. As mentioned in Sec. IV, our second-order 
solution includes the results given in Tomita |Q, Matarrese et al. |T^] and Parry et al. |TB 



although the essential calculation we need in our approach is just the solution of a second- 
order ordinary differential equation by iteration method. Therefore, our approach surpasses 
these others in perfection and simplicity. 

Another advantage of our method remains, which is not mentioned above. Tomita's 
approach is valid only when the absolute value of the density contrast |5| <^ 1 while ours 
does not rely on this assumption, which is the inherent usefulness of the so-called Zel'dovich 
approximation. The gradient expansion technique implies taking the "square root" of the 
metric tensor in order to reproduce the Zel'dovich approximation, while we do not need such 
a trick since we start from the tetrad formalism. 

In our approach the extensions to 7^ 0, A 7^ cases and radiation universes {p = |p) 
are straightforward. These will be the subjects of future investigation. 
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APPENDIX A: GAUGE CONDITION 

The most general gauge transformation to first order is the result of the infinitesimal 
coordinate transformation 
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= + . (Al) 
The changes in the 4-velocity and in the metric tensor are computed from 



which gives to first order 

6gu^ = u^{x^) - u^{x^) = i^u"" - u^^r , (A3) 

If the perturbations are hnear, we can treat scalar perturbations seperate and we write 

e = {T,5'^Lj) (A4) 

for the /c = background, where T = T{x^) and L = L{x^) are arbitrary scalar functions. 
The gauge condition we impose in this paper is the comoving synchronous condition: 

= , c/oo = -1 , ^oi = . (A5) 

These equations must hold for every gauge transformation, so that = ^cfi'oo = ^G%i = 
lead to 

Lj = , T = , T, = . (A6) 
Apart from a trivial constant translation, these are solved to give 

T = 0, L,=L,(x). (A7) 
The change due to the residual gauge freedom is 

6G9ij = -2a'^L^ij{x.) , (A8) 

or if we use Eq. (|3.11| ), 

(5gC,,,(x) = -L,,(x) . (A9) 

Therefore, using the remaining gauge freedom L ., (x), we can choose Cjj(x) = 0. 
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APPENDIX B: COMPLETE SECOND-ORDER SOLUTIONS 



The complete solution for the triad reads 



e 



where ■j/'^^- and 0*^- are the spatially dependent parts of the second-order homogeneous so- 
lutions of Eq. (|2.21| ), (f'^j, v^j and Cj those of the second-order inhomogeneous solutions, 



which come from x \E', x $ and $ x $, and "^^j and O^j are the corresponding 
transverse-traceless parts. 



1. The coupling mode 



We obtain 



k£ 



fB2) 



with 



20 



fin this calculation we use J- = J , which comes from V'\ 

The equation for the transverse-traceless part is 

■' a a 



(B3) 



(B4) 



where 



+ 



Ml 



,3 



2^•^;.<l>•^ + 2$•^fc*•^ - 24^'^;.<l>'^^. - 20^''=$'*j.fc 



i ,k 



(B5) 
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Using the conformal time rj, this is rewritten as 

4 d 



df]'^ ^ T] dr] r] 
The solution is 

^^•(x, 77) = ^ dr' r' ((77 + r')(77 - r') - r'rji^) J dn' r,.(x + x') . (B7) 

2. The decaying mode 

We obtain 

C, = ^(a\5^.-4AV) , (B8) 

where 

A^^. ^ 1 - . (B9) 

The equation for the transverse-traceless part is 



e\ + s-e% - irV^e\ = r^Q\. , (bio) 



where 

Again this is rewritten as 

(9^ A f) 72Q 



Q^2 J 

and we obtain the solution 

729 p-^in I I 



^^■(x, 77) = ^ f^'"'" dr' r' ((477 - r')(77 - r')-' + r'r^-f) j dQ' Q^,.(x + x') . (B13) 
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3. The complete expression of the metric tensor 



The complete metric reads 
20 



9 

1 4 



19*'^,*,,, - 12*'^,*,, + 3 (^(vl/'y - 5. 



+ r 



+ r3 



+ 2^. 



+ (l2^'^^$'^^ - 2^'\^^[^ + 10^''=$'^;;^^) + 9 



k,ij 'rij 



I) 



(B14) 
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